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In this work we investigate the existence of relativistic models for dark matter in the context 
of bimetric gravity, used here to reproduce the modified Newtonian dynamics (MOND) at galactic 
scales. For this purpose we consider two different species of dark matter particles that separately 
couple to the two metrics of bigravity. These two sectors are linked together via an internal (7(1) 
vector field, and some effective composite metric built out of the two metrics. Among possible 
models only certain classes of kinetic and interaction terms are allowed without invoking ghost 
degrees of freedom. Along these lines we explore the number of allowed kinetic terms in the theory 
and point out the presence of ghosts in a previous model. Finally, we propose a promising class of 
ghost-free candidate theories that could provide the MOND phenomenology at galactic scales while 
reproducing the standard cold dark matter (CDM) model at cosmological scales. 

PACS numbers: 95.35.-|-d, 04.50.Kd 


I. INTRODUCTION 

General Relativity (GR) successfully describes the 
gravitational interaction in a wide range of scales and 
regimes, from the solar system size to strong fields in 
binary pulsars and black holes, and most likely will con¬ 
stitute the correct tool for the future gravitational wave 
astronomy [Ij. Up to now, GR has been able to prevail 
against all alternative theories, either scalar-tensor im, 
vector-tensor dQ or tensor-tensor theories, the latter 
comprising massive gravity [ll,[l3, bigravity [13, [3 a-M 
multigravity [l9j | theories. 

In spite of these successes, the extrapolation of GR to a 
broader range of scales — notably, cosmological scales — 
faces important challenges since it relies on the introduc¬ 
tion of a dark sector, composed of dark matter and dark 
energy. The nature of this dark sector constitutes one of 
the most important mystery of contemporary physics. 

The reference model of cosmology today assumes a 
pure cosmological constant A added to the field equations 
of GR to account for the dark energy, and a component of 
non-baryonic dark matter made of non relativistic parti¬ 
cles called cold dark matter (CDM). The best motivated 
candidate for the dark matter particle is the WIMP (2^ . 


‘Electronic address: blanchet®iap.fr 
lElectronic address: laviniah@kth.se 


The model A-CDM is very well tested at cosmological 
scales by the accelerated expansion of the universe, by 
the observed fluctuations of the cosmic microwave back¬ 
ground, and by the distribution of dark matter in large 
scale structures. 

Unfortunately this model does not explain the presence 
of a tiny cosmological constant A. In the prevailing view 
it should be interpreted as a constant energy density of 
the vacuum. However, the unnatural observed value of 
A and the instability against large quantum corrections 
put in doubt its consistency using standard quantum field 
theory techniques. 

Another important concern is that the model A-GDM 
does not account for many observations of dark matter 
at the scale of galaxies, where it faces unexplained tight 
correlations between dark and luminous matter in galaxy 
halos [13, HI] . Primary examples are the baryonic Tully- 
Fisher relation between the asymptotic rotation velocity 
of spiral galaxies and their baryonic mass, and the corre¬ 
lation between the mass discrepancy {i.e. the presence of 
dark matter) and the acceleration scale involved [13, [13| . 
These correlations happen to be very well explained by 
the MOND (Modified Newtonian Dynamics) empirical 
formula [23 - [27l |. The agreement between MOND and all 
observations at galactic scales is remarkable and calls for 
an explanation. On the other hand, MOND has prob¬ 
lems explaining the DM distribution at the larger scale 
of galaxy clusters . 

Many works have been devoted to promoting the 
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MOND formula into a decent relativistic theory. Most 
approaches modify GR with extra fields without invok¬ 
ing dark matter Here we shall be interested 

in another approach, based on a form of dark mat¬ 
ter d la MOND called dipolar dark matter (DDM). 
This approach is motivated by the dielectric analogy of 
MON D 1^ . A first relativistic model was proposed 
in [ 4 ^ |44| and shown to reproduce the model A-CDM 
at cosmological scales. Recently, a more sophisticated 
model has been based on a bimetric extension of GR [1^ 
(see also [1^ for further motivation). In this model two 
species of dark matter particles are coupled respectively 
to the two metrics, and are linked by an internal vector 
field generated by the mass of these particles. The phe¬ 
nomenology of MOND then results from a mechanism of 
gravitational polarization. 

Bimetric theories have been extensively investigated in 
the quest of a consistent massive gravity theory going be¬ 
yond the linear Fierz-Pauli theory. The past decade has 
seen the emergence of a specific theory that avoids 

the appearance of the Boulware-Deser (BD) ghos t to 
any order in perturbations. This dRGT theory [1^, 
has been extended and reformulated as a bimetric theory 
with two dynamical metrics [l3, [13 ■ The theoretical and 
cosmological implications of these theories are extremely 
rich. Notably, cosmological solutions of massive gravity 
theories have drawn much attention (see also the 

references in 0). 

In the present paper we point out that the previous 
model for DDM in a bimetric context [13 , despite the im¬ 
portant phenomenology it is able to reproduce, is plagued 
by ghosts and cannot be considered as a viable theory. 
Nevertheless, this phenomenology (especially at galactic 
scales, i.e. MOND) definitely calls for a more fundamen¬ 
tal theory. We look for a consistent coupling of the dark 
matter helds to bigravity, closely following the restric¬ 
tions made in |52h54| . We thus propose a new model, 
whose dark matter sector is identical to the one in the 
previous model [13, but whose gravitational sector is 
now based on ghost-free massive bigravity theory. As bi¬ 
gravity theory represents essentially a unique consistent 
deformation of GR, we think that the new model will 
represent an important step toward a more fundamen¬ 
tal theory of dark matter d la MOND in galactic scales. 
In a separate paper [53 we work out in more details the 
new model and investigate whether it reproduces also the 
cosmological A-GDM model at large scales. 


II. DIPOLAR DARK MATTER 


A new relativistic model for dipolar dark matter was 
constructed in [13 via a bimetric extension of GR, which 
recovers successfully the phenomenology of MOND. It 
relies on the existence of two species for dark matter that 
couple to two different metrics and an additional internal 


field in form of a vector field, 
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where pbi Pgj Pf are the scalar energy densities of pres¬ 
sureless ordinary matter (baryons) and the two species 
of dark matter respectively, and j'i denote the con¬ 
served currents of the dark matter. On top of the two 
Einstein-Hilbert terms for the g and / metrics, there is 
an additional kinetic term for the effective metric QeS 
and a cosmological constant Aes associated to it (here 
we neglect possible cosmological constants in the g and 
/ sectors). The t/(l) vector field Ag, is introduced to link 
together the two species of dark matter particles and has 
a non-canonical kinetic term yV(A’), with 

T = , (2) 


and Fgv = dgAv — d^Ag. The rich phenomenology and 
physical conse quen ces of this model were studied with 
great detail in [^. For a particular choice of the func¬ 
tion yV it recovers the desirable features of MOND and 
passes the constraints of the solar system. Furthermore 
it agrees with the cosmological model A-CDM at first or¬ 
der cosmological perturbation and is thus consistent with 
the fluctuations of the CMB. 

The effective composite metric GeS was computed per- 
turbatively in [^ and here we show the exact non- 
perturbative solution for this metric. Furthermore, we in¬ 
vestigate the number of gravitational propagating modes 
and the presence of ghost instabilities. The metric C/®® 
was defined in [^ by the implicit relations 


gf. = g^spp^^u = g^sPp.u. ■ 


( 3 ) 


After introducing the matrices G)) = G^gPgp and F))' = 
Gesfpp above relations simply become 


GF = FG = 1, 


( 4 ) 


thus G and F are the inverse of each other. Using this 
fact the form of Geff can be computed. This was done 
perturbatively in [13 with result given by (A8) there. 
Actually the solution of ([3])-(|l]) can be obtained exactly. 
For this, we note that the following relations are true, 

g-^f = {GesG)-^GesF = G-^F = F\ (5) 


This means that we can identify 

F=^/^. (6) 


Thus the exact solution for the effective metric fulfilling 
the relations ([3]) is suggestively 

gf. = 9pp{V^)l- 


( 7 ) 
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Here we will first pay special attention to the conse¬ 
quences coming from the kinetic term and a cosmological 
constant for this effective metric QeS in the action O- 


III. MORE ON EFFECTIVE METRICS 

In the different context of massive bigravity theories, 
interesting proposals for an effective composite metric 
were made in [53 - 1^ . There the form of the effective 
metric was determined by the question of how the cou¬ 
pling of the matter fields to the two metrics of massive 
bigravity behave at the quantum level and whether they 
alter the specific potential interactions of the allowed po¬ 
tential interactions between the two metrics. One par¬ 
ticularly interesting effective composite metric has the 
following form, 

afi, = , (8) 

where a and /3 are arbitrary constants. Defining the 
quantities 9~^fYv as was 

done in (where is shown to be symmetric), it is 
straightforward to see that the determinant of this com¬ 
posite metric corresponds to the allowed potential inter¬ 
actions in massive bigravity, 

det( 5 ®®) = det + ^b"pp) (ag^o- + l3Y^a)] 

= (det g)~^ [det{agfj,^ + ^ 

= (det 5 ) [det(al-I-. (9) 

Thus, the square root of the determinant of 5 ®®, say 
geS = det(g®®), corresponds to 

V-9eS = det(al -b /3X) . (10) 

This is the right form of the acceptable potential inter¬ 
actions between the metrics g and /. Expanding ^—geS 
around a flat background, defining 

/pi/ = (11) 

they correspond to the specific interactions of the form 
4 

y=^=^(a + /3)4-"e„(fc), ( 12 ) 

n—0 

where and the symmetric polynomials 

are defined by (with [• • • ] denoting the trace as usual) 

e^ik) = 1, 
ei(fc) = [k] , 

em = l(W’-3M[)n"]+2[fc»]), 

+8[fc][fc3] -6[fc4]). (13) 


Thus, ^/—geB has exactly the nice structure of the poten¬ 
tial with the special tuning in order to remove the BD 
ghost at any order [ll, [l^ . 

Finally, the relation between the effective composite 
metric QeS proposed in l45ll and the alternative effective 
metric geS proposed in is given by [since tl®® = 
from JT])] 

5p® = a^ffpi/ + 2a/3 a®® + /3Vpi/ ■ (14) 

In other words, QeS = g\/g~^ f does not contain the lin¬ 
ear parts proportional to g and / in ^eff- Unfortunately, 
this will have important consequences as we will see in 
the following section. 


IV. COSMOLOGICAL CONSTANT FOR THE 
EFFECTIVE METRIC 


We will first study the consequences of having in the 
model m the square root of the determinant GeS = 
det(t/®®). Since our non-perturbative solution is t/®® = 

5pp(\/<?“^/)0; the square root of the determinant reads 

V-GeS = \/■ (15) 

Perturbed around a flat background, it corresponds in 
the notation (fTTl) to 

y^Z^=l + l[h + i]+l(^[h + if-2[h^+f]) 

+^{[h + - &[h + i] [h^ +f]+8 [h^ + i ^]) 

+ •••. (16) 


As is immediately seen, does not have the right 

potential structure, in fact it does not even contain the 
right structure for the linear Fierz-Pauli mass term. Any 
Lagrangian that contains this term as a possible poten¬ 
tial interaction between the two metrics has immediately 
the BD ghost at the linear order. Thus the cosmological 
constant for this effective metric or any minimal coupling 
to matter fields via tl®® will reintroduce the dangerous 
ghostly mode. The ghost would come already at a scale 




m^Mp(n7r)^ 

Al 


9 ’ 


(17) 


where Ag = Mpm? and tt denotes the 0-helicity mode. 
This means that the ghost is a very light degree of free¬ 
dom. This immediately kills the possibility of consider¬ 
ing any Lagrangian (independently of all the additional 
terms present in it) that contains y—GeB- 


V. MINI-SUPERSPACE OF THE NEW 
KINETIC TERM 


In the previous section we studied the implications of 
having the cosmological constant for Cl®® and saw that it 
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introduces ghostly interactions between the two metrics. 
In this section we will pay attention to the kinetic term 
\/—QeaTi-eS, where 7?.eff is the Ricci scalar built from C/®®. 
Moreover, we will investigate the allowed number of ki¬ 
netic terms. The first test that such term has to pass is 
the special case of the mini-superspace. The respective 
metrics in the mini-superspace are given by 

dSg = = —n^gdt^ -|- a^da;^ , 

ds^ = fi^udx^dx^ = —njdt'^ + Ojda;^ , (18) 

where Ug, nj and ag, af are functions of the cosmic time 
t only. Consider the following Lagrangian with the three 
kinetic terms 


„ M'^ , _ M? ,- ;-- 

+ -^V^Rf + ^ , 

(19) 

that in the mini-superspace simply becomes 


^eff _ 

'^kin 




SMgjjaeffdgff 


n/ 


^eff 


( 20 ) 


Following the prescription (|3]) we obtain Ugg = y/rigUf 
and Oeff = y/agQf. We compute the conjugate momenta 
for the scale factors and get 


Pa = -dM^alHg - [^^{HgUg + HfUf) , 

z y rigUf 

Pf = -GM]a}Hf - ^M^gag [^^{HgUg + HfUf) , 
z y rigrif 

( 21 ) 


interactions. Because of their very different structures 
there is no hope for cancellations between these terms. 

Taking the limit when Mf —^ 0 of the Hamiltonian 
(1^ results in 


n 


eff 

kin 




1 

1^ 



-PfO-ffng 

M! 


+ 


Ap)ag^Ja^^/n^ \ 




(24) 


As one can see, even in this limit the Hamiltonian is 
not linear in the lapses, so that the variation of the 
Hamiltonian with respect to the lapses gives rise to equa¬ 
tions of motion that depend on the lapses and hence the 
constraint equation is lost. Therefore, the kinetic term 
V~Ges'R'eS introduces ghostly interactions already in the 
mini-superspace independently of the number of present 
kinetic terms. 

An interesting question to address at this stage is 
whether or not the mini-superspace can be made ghost- 
free by considering the kinetic term for ^eff that was pro¬ 
posed in [s^. Since the determinant of geS corresponds 
to the right ghost-free potential interactions between two 
metrics, the kinetic term for ^eff might behave better than 
that for t/eff. Thus, consider as next the Lagrangian with 
the alternative three kinetic terms 

Ain = —V^Rg + —V^Rf + —^V-gesReS, 

(25) 

where i?eff is now the Ricci scalar of the metric g®®. In 
the mini-superspace this becomes 


where Hg = and similarly Hj are the conformal 

Hubble factors. Now we can perform the Legendre trans¬ 
formation to obtain the following Hamiltonian: 

■^kfn = ^{alnf(^M^sply/d:^ng + AM^p^af 

+ ‘^Pg^gafUg (-M^sPf^agafUf + 2Mjpgafy/n^'^ 

+ V ^g^f^g^f I ! (22) 

where we defined the shortcut notation for convenience 


^eff _ 

■^kin 


ZM^agh 


2 

9 


iMjafa) 3M^s'a,sals ^ ^26) 

n/ heff 


with this time heff = aug + jdnf and heff = aag + f3af. 
The conjugate momenta for the scale factors are now 

Pg = -6(^M^alHg + ^^^^^^{aagHgUg + PafHfUf)^ , 

Pf = -d(M]a)Hf + ^M^^I^aagHgUg + PajEfUf)^ . 

(27) 


2 = -l2(^MlgM]a),/a^ng + M'^sMgal./a^nf 

^AM]Myga)^n^) . (23) 

The Hamiltonian is highly non-linear in the lapses ng and 
nf. Since there is no shift over which we have to inte¬ 
grate, this is an immediate sign that these three kinetic 
terms have the BD ghost degree of freedom already in 
the mini-superspace (see [s^ for an introduction to con¬ 
strained hamiltonian systems). Thus, one has to avoid 
the two very bad contributions in form of (i) the cosmo¬ 
logical constant term for GeS, and (ii) the kinetic term 
VSesR-eS — both these terms correspond to ghostly 


Thus, the Hamiltonian is given by 

^kfn = ^{agnf[-a{M^pj + M^g{apf - ^Pgf)ng 
- M^pjj3nf] + afUg [-M^fgUng - P{M]fg 
+ M^do^pf - I3pgf)nf] I , (28) 

where Q stands for 

2 = l2(^Mjaaf[{Mg + M^ga'^)ag + MggaPaf]ng 

+Mgl3ag [Mfa / -b Aiydes] Uf^ . (29) 
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Again, the Hamiltonian is highly non-linear in the lapses. 
The problem comes from the fact that we have too many 
kinetic terms. Indeed we see immediately that the only 
way of having linear dependence in the lapses in the mini¬ 
superspace (and hence getting rid of the BD ghost) is if 
we take either the /3 —>■ 0 limit — this would simply 
correspond to having only the standard kinetic terms for 
the g and / metrics — or the M/ —> 0 limit, 

■yeS I ^ ^esPjo^Pf - PPgfafng 

a{Mlp) + M^giapf - PPgf)ng + M^ppuf 

As we see, the Hamiltonian becomes linear in the lapses 
when we remove for instance the kinetic term for the / 
metric. Thus, the only way of having a healthy mini¬ 
superspace is if we restrict the kinetic Lagrangian to be 
either Mgy/—gRg + Mjy'—fRf which are the standard 
ghost-free kinetic terms, or M^^/^Rg + M^g^y—gesReS- 
In a symmetric manner we could also remove the kinetic 
term for g and hence Mj-y—fRf+M^g^y—gegRes would 
be also perfectly valid. In summary, one should restrict 
the theory to have not more than two kinetic terms in 
order not to reintroduce the BD ghost. 


VI. DIPOLAR DARK MATTER IN 
GHOST-FREE BIMETRIC THEORY 


The dark matter model proposed in is therefore 
non-viable, but nevertheless points toward an interesting 
connection between dark matter at small galactic scales 
(interpreted as DDM) and bimetric gravity. Based on 
our previous analysis, we would like now to propose the 
following new model for dipolar dark matter based on 
ghost-free bimetric theory, 


-fK 


— \/ g ^ 2 ^ ^9 Pi' Pgj ” 1 " 7 2 


-fM'i 


“b "n/ l?eff 


+ Ag{j^g - fj) +W{X) 


( 31 ) 


where the ghost-free potential interactions are defined 
by the metric ([U [they take the form (fT^ -dT ^ when 
expanded around a flat background], and where the ki¬ 
netic term of the vector field is now constructed with the 
metric g®®, 

^ ■ (32) 

As was shown in , the matter fields can separately 

couple to either the g metric or / metric without invoking 
the BD ghost. Additionally the matter fields can couple 
to the effective composite metric geS which is ghost-free 
in the mini-superspace and in the decoupling limit. 

Here we propose to couple the ordinary baryonic fields 
with mass density pb to the standard g metric while cou¬ 
pling the two species of dark matter with densities Pg 


and pf separately to the g and / metrics respectively. 
Furthermore, in order to link together the two species 
of dark matter particles, we consider a vector field Af_i 
that minimally couples to the effective metric peff- The 
vector field plays the role of a “graviphoton” since it is 
generated by the mass currents and of the parti¬ 
cles. The presence of this internal field is necessary to 
stabilize the dipolar medium and is expected to yield the 
wanted mechanism of polarization. 

The model (l5T]l fulfils the restrictions coming from our 
previous analysis, as it contains no more than two kinetic 
terms (in particular the problematic kinetic term for 
is absent), and the potential interactions between the two 
metrics coming from the square root of the determinant 
of peff correspond to the ghost-free prescriptions. 

However, one needs also to be careful with the assump¬ 
tions on the dark matter helds and their currents. Indeed 
the ghost could still be present in the matter sector.^ 
Since the dark matter fluids that live on the g and / 
metrics directly couple to Ag that lives on the peff metric, 
there is a priori the danger of having the ghost present 
due to the interaction term in the matter sector. We will 
investigate this question with great detail in and see 
if for a specific choice of the dark matter fields the ghost 
can be maintained absent. For this we shall perform the 
decoupling limit analysis of our new model (I3ip includ¬ 
ing the matter sector and study the required amount of 
initial conditions. 

Finally, the model (iMl) should share the nice proper¬ 
ties and phenomenology of the model proposed in 
and therefore provides a promising road for a relativis¬ 
tic dipolar dark matter model to be investigated. The 
MOND phenomenology, the PPN parameters and the 
cosmol ogy of this model will be studied in a separate 
paper |55l |. 


VH. CONCLUSIONS 

We explored the possible candidates for relativistic 
dark matter models in bimetric extensions of General 
Relativity, that hopefully will provide modified Newto¬ 
nian dynamics (MOND) at galactic scales while giving 
rise to an expansion at cosmological scales. A promising 
road comes from the ghost-free constructions of dRGT 
massive gravity [l5l llbl j where the interactions between 
two metrics are tuned in a way that the Boulware-Deser 
ghost remains absent. Furthermore, the important stud¬ 
ies of possible consistent couplings to matter fields Hi- 
[s^l are beneficial to us, since for the model to work, we 
have to consider two different species of dark matter par¬ 
ticles that couple separately to the two metrics while an 
additional internal vector field couples minimally to an 
effective metric built out of the two. The vector field 


i We are grateful to Claudia de Rham for discussions on this. 
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links together the two sectors of the dark matter parti¬ 
cles and plays a crucial role for gravitational polarization 
and MOND [ii,|4i]. 

For the ghost absence the question of allowed kinetic 
interactions is mandatory. We showed that the kinetic 
Lagrangian containing three kinetic terms immediately 
gives rise to the introduction of the ghost and we there¬ 
fore concluded that only two kinetic terms are allowed. 

In a future work [55| . we will study in detail the co¬ 
variant equations of motion of the new model, derive the 
non-relativistic limit and see if the polarization mecha¬ 
nism for dark matter works in the same way as in the 
originally proposed model. We will investigate in detail 
the possible danger of ghostly interactions in the matter 
sector and constrain further the model. We intend also 


to check if the parametrized post-Newtonian parameters 
are close to the ones of GR in the solar system, and to 
investigate the cosmological solutions in first order per¬ 
turbations. 
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